Quantum state preparation and macroscopic entanglement in gravitational-wave 

detectors 
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Long-baseline laser-interferometer gravitational-wave detectors are operating at a factor of ~ 10 
(in amplitude) above the standard quantum limit (SQL) within a broad frequency band (in the sense 
that A/ ~ /). Such a low classical noise budget has already allowed the creation of a controlled 
2.7 kg macroscopic oscillator with an effective eigenfrequency of 150 Hz and an occupation number of 
~ 200. This result, along with the prospect for further improvements, heralds the new possibility of 
experimentally probing macroscopic quantum mechanics (MQM) — quantum mechanical behavior 
of objects in the realm of everyday experience — using gravitational-wave detectors. In this paper, 
we provide the mathematical foundation for the first step of a MQM experiment: the preparation of 
a macroscopic test mass into a nearly minimum-Heisenberg-limited Gaussian quantum state, which 
is possible if the interferometer's classical noise beats the SQL in a broad frequency band. Our 
formalism, based on Wiener filtering, allows a straightforward conversion from the classical noise 
budget of a laser interferometer, in terms of noise spectra, into the strategy for quantum state 
preparation, and the quality of the prepared state. Using this formalism, we consider how Gaussian 
entanglement can be built among two macroscopic test masses, and the performance of the planned 
Advanced LIGO interferometers in quantum-state preparation. 

PACS numbers: 42.50.Xa, 42.50.Lc, 03.65.Ta, 03.67.Mn, 04.80.Nn, 95.55.Ym 



I. INTRODUCTION 

An international array of first-generation long-baseline 
laser interferometric gravitational-wave (GW) detectors 
(LIGO [H, VIRGO H, GEO [| and TAMA (3) are 
reaching their design sensitivities. These Michelson in- 
terferometers have been built to measure GW-drivcn rel- 
ative length changes (within a detection band from 10 Hz 
to 10 kHz) between the mirror-endowed test masses which 
are hung as pendulums with an eigenfrequency far below 
the detection band. Resonant cavities are used to en- 
hance the sensitivity by increasing the circulating optical 
power and the signal storage time. In Michelson inter- 
ferometers usually the total differential mode of motion 
between the arm cavity mirrors, in the following always 
represented by the position operator 

X = (xitMe ~ ^ETMe) ~ (^ITMn ~ ^^ETMn) , (1) 

(for the nomenclature see Fig.[T|) is measured via a ho- 
modyne detection of the modulation fields (also called 
side-band fields) leaking out at the dark port of the in- 
terferometer. Current GW interferometers are already 
quantum limited at high frequencies by the shot noise. 
Next generation interferometers, such as the planned Ad- 
vanced LIGO detector [1] (cL Fig.[T]), are expected to be 
quantum noise limited at nearly all frequencies in the 
detection band. Advanced LIGO will therefore operate 
at its free mass standard quantum limit (SQL) ^ lZi, iS] 
at which the back-action noise is comparable to the shot 
noise level. The position-referred spectral density of the 



SQL at the sideband frequency O is given by 



^SQL(r!) = 



2h 

m 



(2) 



The SQL is the minimum noise spectrum achievable by 
a linear quantum measurement of the position with un- 
correlated sensing and back-action noise. Here m is the 
reduced mass of all arm cavity mirrors, or 1/4 of the 
individual mirror mass. 

Improvement of sensitivities to GWs beyond Advanced 
LIGO will require to surpass the SQL significantly in a 
broad frequency band. Various conceptual strategies ex- 
ist for building interferometers with broadband sub-SQL 
quantum noisei, [H, [Til, [H, El, E3, El, El, E3, El, [ll, 
[20l . [2ll . [H, [H, [U , while there is also much effort devoted 
to lowering classical noise below the SQL. For exam- 
ple, (i) the CLIO interferometer is currently being cooled 
down to a few tens of Kelvin, and has a theoretical ther- 
mal noise budget below the SQL [ll], (ii) non-spherical 
mirrors are being developed that support non-Gaussian 
modes which average better over mirror surfaces and are 
thus less susceptible to thermal noise [H, [13, [11], (iii) 
coating structures are also being optimized for coating 
thermal noise [H, [sl] . 

This paper, however, is devoted not to the improve- 
ment of the detector's sensitivity to GWs, but to the 
study of quantum mechanical behavior of its test masses. 
Recently, it has been reported by the LIGO Scientific 
Collaboration in Ref. [3l| that a certain (controlled) 
mode of the mirror's differential motion in the LIGO 
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FIG. 1: Schematic plot of the planned Advanced LIGO 0] 
interferometer: a power- and signal-recycled Michelson inter- 
ferometer with cavities in the arms and a homodyne detection 
scheme at the dark port. 



detector located in Hanford was experimentally cooled 
down to lAfiK, with an cfFcctive occupation number of 
around 200. Furthermore, there is a number of other 
experiments conside ring different smaller-scale mechan- 
ical structures [H, [sj [13, [M HI, [ll. The goal of 
such cold damping [38[ experiments is to reach the os- 
cillator's ground state. As a real pure quantum state is 
approached, the semiclassical model as used in the above 
references will certainly break down and the quantum 
noise effects in the measurement process have to be in- 
cluded. 

One aim of this paper is at providing a mathemat- 
ical foundation guiding future experimental efforts of 
quantum-state preparation — one that is straightforward 
to apply to experimental situations. Note that our gen- 
eral formulation applies not only to gravitational-wave 
detectors but to the whole bunch of experiments related 
to quantum state preparation. Quantum mechanically, 
we consider a problem in which an object is being con- 
tinuously measured by the optical field, while it is si- 
multaneously subject to noisy forces. For such problems, 
the stochastic master equation (SME) [3§| is a readily 
available tool to simulate quantum-state preparation in 
quantum mechanical systems; a Riccati equation associ- 
ated with the SME provides the remaining uncertainty 
of the object, when all measurement data are taken into 
account. Formally, this approach only treats Marko- 
vian systems (since only Markovian systems allow us to 
project the optical field being measured and trace off the 
noise fields at every infinitesimal time step, resulting in 
a closed-form evolution of the test-mass density opera- 
tor), while non-Mar kovianity is prevalent in experiments 



such as in GW detectors: virtually none of the noise 
sources are white, and we have the added complexity 
that the cavity mode often couples strongly with the test 
masses, i.e. it cannot be adiabatically eliminated (soj . and 
must be evolved together with the test masses. Since we 
only consider linear systems with Gaussian noise, and 
we only care about the test masses' state after the ini- 
tial transient has died down, the SME and the Riccati 
equation, which also characterizes the exact way of the 
transient decays, are not entirely necessary. Instead, we 
have found that a Wiener filtering approach, in which 
the measurement data is filtered with the optimal causal 
Wiener filter to obtain instantaneous optimal estimates 
for position and momentum of the test masses [i^l , suf- 
fices, and is most straightforwardly connected with ex- 
perimental calibration of the system. An example of the 
power of the Wiener filtering approach has already been 
demonstrated in Ref. [4l| . where an optimal controller 
that yields a steady quantum state with minimum un- 
certainty (or von Neumann entropy) has been derived for 
a general linear (Markovian as well as non-Markovian) 
quantum measurement process. 

A direct application of the Wiener filtering approach is 
to explore how quantum a macroscopic test mass can be 
prepared in a certain environment, which is the second 
aim of this paper. We measure the purity of the quantum 
state of an individual test mass (or a single mechanical 
degree of freedom) through the Heisenberg uncertainty 
product or the following quantity 



-I 



^JVxx Vpp V^p , 



(3) 



which is unity for a pure state. This quantity can also be 
converted into an effective occupation number. Here Vxx , 
Vpp and Vxp are the second-order moments of position 
and momentum of the Gaussian state. We will show that 
a low classical noise budget which is completely below the 
SQL for a broad frequency band allows the quantity U 
to become close to unity. In particular for simple sys- 
tems with a total classical noise spectrum S'ci(i^) which 
is dominated by a white sensing noise and a white force, 
we have obtained the simple relation 



U + mm i — }■ 



(4) 



But we will also explore how a realistic noise budget for 
the planned Advanced LIGO detector — as well as an ex- 
tension of Advanced LIGO with plausible improvement 
— can best be taken advantage of through an optimized 
optical configuration that minimizes U. When two inde- 
pendent mechanical degrees of freedom are considered, 
the formalism, which we present in this paper, has al- 
ready been applied to show that the production of quan- 
tum entanglement between the macroscopic end mirrors 
is possible for sub-SQL laser interferometers [i^ . 

An experiment testing MQM should be divided into 
different stages which are separated in time: a prepa- 
ration stage, where the test mass will be continuously 
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observed; an optional free-evolution stage; and a veri- 
fication stage. One will need to collect statistics from 
a huge number of identical trials as it is required from 
quantum mechanics. This present paper is the first one 
of a paper-series: this one deals with the preparation of 
macro scop ic conditional quantum states; while a second 
paper [42| will study the verification of such macroscopic 
quantum states. 

This paper is organized as follows: in Sec.|TT] we will 
briefly review the theoretical basics of Wiener filtering. 
In Sec. lIIII we will study analytically the conditional vari- 
ances of the simplified model using only Markovian dy- 
namics. We will introduce a flexible homodyne detection 
angle and input squeezing. In Sec. lIVI we will extend the 
analysis done in Ref. [40| about macroscopic entangle- 
ment. In Sec.|V] we will study test-masses in a cavity 
with finite bandwidth and we will treat more realistic, 
non-Markovian noise sources. These preliminary studies 
result in an investigation of quantum-state preparation in 
Advanced LIGO in SecEB Finally, in SecEII] we will 
summarize our main conclusions. 



II. WIENER FILTERING 

For systems under continuous measurement, the con- 
ventional approach is to describe the joint system- 
measurement-data evolution using a stochastic master 
equation (SME) [H, H^, [i^, |3l , which is a set of stochas- 
tic differential equations that simulates the joint evolu- 
tion of the system's conditional density matrix and 
measurement data y{t). As the simplest example, for a 
harmonic oscillator with position x being measured con- 
tinuously by a Markovian measurement device (one that 
has uncorrelated measurement noise at different times 
and constant measurement strength), the SME reads. 



+ {{x, ~ 2 {x) dW , (5) 



dy:^-{x) dt + dW/V2, 



(6) 



with a coupling constant a. Here the conditional quan- 
tum state 4^ is defined as the projection of the joint 
system-device quantum state to the sub-space in which 
the readout observable y has definite values of y{t') = 
y{t') y < t' < t. The Wiener increment dW describes a 
stochastic process that simultaneously drives the condi- 
tional quantum state and the measurement data, where 
both are stochastic processes. Different realizations of 
dW correspond to different possible scenarios that could 
take place in reality. 

In practice, it is not enough to only describe the 
stochastic process, we need to be able to obtain the con- 
ditional quantum state at any given time i, based on the 
system's initial quantum state p'^°"'^(0) and measurement 
results {y{t'), < t' < t}. This corresponds to a, filtering 



problem in classical stochastic calculus. The probability 
distribution of any state variable x is simply the condi- 
tional probability 



P[xit)\{yit'),0<t' <t}], 



(7) 



while the conditional expectation of x can be written as 
a functional over {y{t'), < t' < t}, 



^ E[x{t) I {y{t'), 0<t' < t}] 



(8) 



For linear systems with Gaussian states only expecta- 
tions of quantities linear and quadratic in the variables 
are needed. The former can be obtained through a linear 
filter over y, while the latter can be obtained by solving a 
time-domain Riccati equation. If the measurement pro- 
cess has started for sufficiently long time — much longer 
than the time constant of transients — then the filters 
over y as well as the second-order moments are station- 
ary. Theycan be obtained through the theory of Wiener 
filtering [4y| . This will be the situation that we will con- 
sider in this paper. 

Although the most general theory of quantum filtering 
is quite distinct from classical filtering [47|, we will sim- 
plify the situation through two steps, and return to an 
essentially classical filtering problem. First, let us con- 
struct a model of the quantum measurement process. Let 
y(t) be the Heisenberg operator of the measurement out- 
put (e.g., a particular quadrature of the out-going optical 
field) and s{t) any system observable. Then the principle 
of simultaneous measurability, i.e. the observable y(t') at 
different times can be measured individually to arbitrary 
accuracy without imposing any fundamental limits, and 
the principle of causality, i.e. measurement observable at 
present does not respond to future forces to be exerted 
onto the system, will dictate that [i^ 



0=[6(t), y(0] yt>t'. 



(9) 
(10) 



Note here that o could be any any state variable such as 
momentum and position, or even the density operator of 
the system being measured. This means that any filtering 
of the operator y{t') for < i' < i can be considered as a 
classical process, i.e. treated with classical linear control 
theory, as long as the state of the system at i, or p{t) is 
considered. Here we write formally that 



'-cond 

P{v(t'),0<t'<t} 



'P[v{t')=vit').t'<t] P V[y(t')=y(t'),t'<t] 

[Pm')=y(f),f<t\ P^[y(t')=y(t').t'<t]] 



,(11) 



where V projects onto the subspace on which the mea- 
surement operator takes the measured value. Henceforth 
the dependence of ff°^^'^ on {y{t'),t' < t} will not be 
written explicitly, as it is done in most of the literature. 

Second, we restrict ourselves to stable linear systems 
with Gaussian states and to later times when the initial 
states are no longer important. In this case only first- 
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and second-order moments in linear observables need to 
be considered. We also assume all linear observables to 
have zero unconditional expectation values. Suppose we 
have observables xi with I = 1, 2, ... n, and proceed as in 
classical Wiener filtering. The key is to decompose the 
operator xi into 



xiit) 



dt' Ki{t~t')y{t')+Ri{t) 



with 



{Ri{t)y{t')) ^ir pRi{t)y{t') 



0. 



(12) 



(13) 



If this can indeed be done, it can then be shown that 



^^"^(0 =tr [r^'^ut)] = 

and 



v; 



— ' ^ {xi{t)S^m{t))r^ 



I sym 



sym 



dt' Ki{t-t')y{t') (14) 

{xi{t)r°^^{x„,{t)r°^^ 

(15) 



In this way, the conditional expectations of the linear 
variables are given as linear functionals of past measure- 
ment data {y{t'), t' < t}, and the conditional variances 
as steady state constants. Please see Appendix|A] for a 
more rigorous justification. 

We will now try to obtain Ki and VJ™'^'* in terms of un- 
conditional correlation functions, or cross spectra among 
system observables xi and the output y. Eq. ([13]) leads 
to the Wiener-Hopf equation 

C,,v{t - t") - f dt' Klit - t') Cyy{t' - t") = 

J — oo 

yt" < t . 

(16) 



Here 



Cab{t-t') = {a{t)kt'))sym 

a{t)b{t') + b{t')a{t) 



tr 



(17) 



stands for the symmetrized time-domain two-point cor- 
relation function between two arbitrary Heisenberg oper- 
ators a{t) and h{t). 

If we suppose that Ki{t) = for t < 0, i.e. we make 
sure that the filter is a causal function, we can rewrite 
Eq. IHl) as 



C,,y(t)- dt' Kl{t')Cyy{t-t')=Q > . (l8) 



In Fourier domain the condition in Eq. is satisfied if 
the function 



L{n) ^ s^,y{n) ~ Ki{n) Syy{n) 



(19) 



is analytic in the lower-half complex plane while the 
Fourier transform of the filter function Ki{Q) is analytic 
in the upper-half complex plane. Furthermore, L(f2) has 
to vanish at infinity, because Eq. (|18p must also be valid 
at t = 0. These conditions uniquely define Ki{Q). We 
have denoted with Sab{^) the single-sided (cross-) spec- 
tral density among the operators a and 6, related to the 
two-point correlation as 



Cab{t) 



2^ 



-iVtt 



(20) 



The Wiener-Hopf method provides the solution for Ki{^) 
as 



1 



4(17) 



Sy (Q) 



(21) 



where we have split Syy{Q) ~ Syi^) Sy {^) in such a way 
that (s~) and its inverse are analytic functions in the 
upper-half (lower-half) complex plane. Because Syy{fl) 
is in general a rational function of fl^ with real coeffi- 
cients, we expect that Sy{fl) = [Sy{fl*)]*. In addition, 
[F(ri) . . .]+ stands for taking the component of a func- 
tion whose inverse Fourier transform has support only 
in positive times. Operationally, this could be obtained 
either by decomposing F{Q) into 



Fm = E 



(22) 



and only keeping terms whose Clk has negative imaginary 
parts, or by inverse Fourier transforming F{fl) into the 
time domain, eliminate the positive-time component, and 
then Fourier transform back. [Note that both approaches 
will become ambiguous when F{Q) does not approach 
zero when +oo.] Inserting Eq. ((2T|) into Eqs. ([12]) 

and (1151). we obtain 














+ 


. Sy . 





(23) 



A more transparent understanding of the filtering 
problem can be obtained when we formally define the 
causally whitened output as 



1 



4 m 



(24) 



Because 1/SjJ' is analytic in the upper-half complex plane, 
z can be written as an integral over the history of y. The 
random process z has a white spectrum. Moreover, we 
can write 



In other words, 



dt'C,,,it-t')zit') 



(25) 



(26) 
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and 



cond 



Vlr^ 



dt'C^,,{~t')C^^4^t'). (27) 



V, 



Clearly, Eqs. (|26p and ((27)) are simply continuous ver- 
sions of linear regression over a set of independent ran- 
dom variables. 

In this paper with tools of Wiener filtering, we no 
longer need to write down SMEs. In fact, for realis- 
tic systems with multiple colored noise sources and non- 
Markovian dynamics, SMEs can only be obtained by in- 
creasing the dimension of the problem, which will defi- 
nitely become cumbersome. 



III. MEASUREMENT CONFIGURATIONS 
WITH MARKOVIAN NOISE 

A. General discussion 

Let us start our discussion generally with an ab- 
stract continuous linear Markovian measurement process, 
which monitors the center of mass position of a simple 
harmonic oscillator. The Heisenberg equations of motion 
in the frequency domain can be written as 



y = Z + x, 
X = R^^[n) F, 



(28) 
(29) 



where the two noise operators Z and F both have a white 
spectrum. In the time domain, Z and F are white noise, 
with two-time correlation functions proportional to delta 
function. Such statistical characteristics make the mea- 
surement process a Markovian one. Thus, we assume 
that Z and F have the single-sided (cross-) spectral den- 
sities Szz > 0, S'ff > and Szf G which satisfy the 
Heisenberg relation of the measurement process Q 



Szz Sff 



o2 



ZF 



M > 1 



(30) 



which arises from the requirement that the Heisenberg 
operators of the out-going field at different times must 
commute, and guarantees that the level of back action is 
just enough to enforce the Heisenberg uncertainty rela- 
tion of the test mass . In the case of Gaussian noise 
only, we have = 1 if and only if the measurement pro- 
cess is purely quantum. 

The linear response function of a damped harmonic 
oscillator is given by 



1 



TO {^'^ + i^m^ - t^m) 



(31) 



with the eigenfrequency ujm, the damping rate 7„ ^ 
and the mass to. Then we can easily assemble the single- 



sided spectral densities of the measurement process as 

Syyin) = m^\R,,in)\^ {n-^ ~ n^{2qi - 7^) + qi) Szz , 

(32) 

s^yin) = -m\RU^)\^ (n" - {-Iran - - ^) Szf , 

V JZF / 

(33) 

sun)^\R,,{n)\^SFF, (34) 



where we have defined the coefficients 



_ 2 , Szf 



mSzz ' 



mbzz m^Szz 



(35) 
(36) 



Both, 51^2, have the dimension of frequency squared, 
while q2 is always positive. We also note that |gi| < q2. 

From Eq. ([32]) we can recover a quantum limit of the 
measurement process: if we have Szf = 0, the spectral 
density of the measurement noise satisfies 



Syy{n) >2^SzzSFF\R.xm\ > 2^h\R,,{n)\ 

= VM^SQL(f^).(37) 



For a free mass we can then recover Eq. ([2]) from Eq. (|37p . 
Note that for the free mass the first inequality sign in 
Eq. (P7|) becomes an equality sign at VI = Vlq, where the 
measurement frequency is defined by 



= 



Sff 



izz 



(38) 



Therefore, the measurement frequency fig is the fre- 
quency at which the noise spectral density of a Marko- 
vian quantum measurement process (with Spz = 0) ap- 
proaches most its free mass SQL. 

Now it is straightforward to derive the conditional vari- 
ances assuming p = —imSlx and using the Wiener fil- 
tering method, i.e. Eq. ([23]) . For this it has actually been 
crucial to first spectral factorize Syy{il.). In the perfect 
oscillator limit with 7™ — *■ the conditional covariance 
matrix can be put into the following form 



v = vmD 



f I 2(12 jq2 - qi \ 
J 



91+92 V 92 + 91 



92 - 91 



2 



92 



\ V 92 + 91 V 91 + 92 
where we have defined 

/ 

D = 



(39) 




(40) 
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Eq. (|39p gives us the most general covariance matrix of 
the conditional Gaussian state of a lossless harmonic os- 
cillator under any linear Markovian position measure- 
ment. In the covariance matrix V, the matrix D sets 
the scale of the quantum fluctuations to be comparable 
to those of the vacuum state of a harmonic oscillator with 
eigenfrequency ^Jq2- Depending on the ratio (71/(72, which 
always lies between —1 and 4-1, the noise ellipse of the 
vacuum state is deformed into one where position and 
momentum are correlated (unless if (71/^2 ~ 1) while the 
area is conserved. This corresponds to a unitary trans- 
formation among Gaussian states. Then p. finally en- 
larges the noise ellipse with a uniform factor, converting 
the pure state into a mixed state (unless if /i = 1). We 
have found that the conditional variances completely co- 
incide with those obtained from SMEs. With Eq. (|39p we 
can for instance easily reproduce Eqs. (2.8a)-(2.8c) from 
Ref. [13. 

A Gaussian state is pure if and only if its uncertainty 
product is Heisenberg-limited. Therefore, it makes sense 
to quantify the purity of the conditional state by its un- 
certainty product, here given by 




\ - - - V with Hq = Hp / 3 

Vs5' with Hq = VHf ily, 

VS? with Ha = 3 Hx 

\ 

\ 

.V - - - 

\ 

s 

\ 

\ 



n [Hz] 

FIG. 2: Example noise spectral densities (in arbitrary units) 
of a Markovian measurement process observing a free mass: 
the quantum noise spectral density at different values of the 
measurement frequency fi, as well as the spectral densities 
corresponding to a classical force noise (gray line) and a sens- 
ing noise (gray line) are marked in the plot. The gray shad- 
owed region marks the classical noise SQL beating which we 
have chosen to have ili/r^i? = 5. 



det V = 14a: K, 



V 



xp 



(41) 



which is identical to the purity of the measurement pro- 
cess. This simply shows that in the Markovian case, any 
measurement will produce a pure conditional state of a 
lossless harmonic oscillator if and only if it is a quantum 
measurement. Moreover, the uncertainty product is even 
independent of the system's mechanical properties such 
as the oscillator's mass and eigenfrequency. Note that in 
the Appendix|B]we have introduced how the uncertainty 
product is related to an effective occupation number AfeS- 
It turns out, when using Eq. (|B3[) . that here the effec- 
tive eigenfrequency for the effective occupation number 
is given by Woff = ^92- 

The covariance matrix in Eq. (j39[) becomes obviously 
diagonal and the correlation between x and p in the con- 
ditional state vanishes if and only if qi = (72. But this 
is strictly forbidden due to the Heisenberg uncertainty 
principle. However, in a certain limit they can become 
very close [4l| . With a higher difference in qi and (72 the 
correlation in x and p increases — and in turn also the 
uncertainty product. 

In order to obtain the conditional state, as given in 
Eq. (|39p. the measurement data has to be filtered in real 
time using the Wiener filter functions for position and 
momentum which are given in the frequency domain by 



qi 



(42) 



K,in) = im ((72 - ^ ' (43) 



where fli^2 = i±VQ^ 



^qi - W12 - 9i)/V2 and fls 
qi. Note that the poles of the 



i/%/2(. 

Wiener filter fli^2 are actually equal to the zeros of the 



measurement's output spectrum Syy{fl), which in turn 
correspond to the frequencies of maximal sensitivity and 
are therefore easy to find. 

Let us now have a closer look on the noise model we 
will use throughout this section. We will only consider 
quantum measurement processes and categorize the noise 
into two groups: (i) the one which is a result of the mea- 
surement process itself will be denoted by quantum noise; 
and (ii) the additional noise will be called classical noise 
which does not directly arise from the measurement pro- 
cess, usually has no correlation in Z and F does not have 
to satisfy Eq. (|30p . Then the two noise sources combine 
as Sff ^ Sl,p + S'j^p and Szz = •S'l^ + ^zz- 

The quantum noise is dominated at high frequencies 
by shot noise which is covered by and at low fre- 
quencies by back-action noise which is covered by Spp. 
The latter one is represented by the radiation-pressure 
noise in the case of a measurement with light. If both 
are uncorrelated, i.e. S'^p = 0, they result in the SQL. 
Then the quantum noise spectral density is limited from 
below by the free-mass SQL as shown in Eq. ([37|) and also 
in Fig. [2] The quantum noise touches the free-mass SQL 
at the frequency = Qq. 

The other noise source can also be divided into two 
parts: a classical force noise Spp is added to Sff which 
acts directly on the center of mass of the measured ob- 
ject, and is in real interferometric experiments due to 
for instance seismic noise or thermal noise in the suspen- 
sion of the mirrors. The classical sensing noise S"^^ is 
only a pseudo motion of the measured object and may 
be due to the following reasons: (i) on. the one hand due 
to thermal fluctuations of the mirror's shape as for ex- 
ample mirror internal thermal noise which makes only 
the mirror surface move with respect to its center of 



7 



mass; or (ii) on the other hand be due to optical losses; 
or (in) due to photo-detection inefficiency. Therefore, 
our sensing noise is somehow generalized from what is 
conventionally understood when using the term "sensing 
noise". Note that for detecting gravitational waves, only 
the total noise matters, yet for studying quantum-state 
preparation, it is important to make distinctions between 
sensing and force noise, and between quantum and clas- 
sical noise. Throughout this section we will assume both 
classical noise sources to have a white spectrum. Then 
they can be characterized by the frequencies, Qp for the 
force noise and Qx for the sensing noise, at which their 
noise spectral density intersects with the free-mass SQL 
(cf. Fig. [21). These classical noise frequencies are defined 
by the following relations 



S' 



FF 



2hmn 
2h 



F ; 



zz 



(44) 
(45) 



Note that the classical force noise together with the clas- 
sical sensing noise can open a window in which both 
are below the free-mass SQL as indicated by the gray- 
shadowed region in Fig. [21 The total classical noise is 
then even completely below the SQL in a certain fre- 
quency band, if and only if the classical noise sources 
satisfy 



ccl ccl ^ fc2 

^FF '^zz ^ " 



<j4> nx/nF>2, 



(46) 



which turns into a constraint for the classical noise fre- 
quency ratio. If Eq. (|^^ holds, the classical noise is equal 
to a minimum of (2rii?/f2.^) times the free-mass SQL at 
the frequency Q = fid = which reads 



5ci(r>) I _ ^ci(»ci) _ 2^!^ 

^ \SsQl{^)} <S'sQL(^^cl) 



(47) 



Since here the classical noise has the largest separation 
to the SQL, we can understand the factor {2flp/nx) as 
the classical-noise-SQL-beating factor. 



B. Very low finesse cavity and free mass scenario 
witli vacuum input 

Let us consider the simple situation of a laser beam 
incident on a suspended mirror, where the output field 
is monitored by a perfect balanced homodyne detection 
at a frequency-independent angle This corresponds 
to a mirror in a cavity with infinitely large bandwidth 
(or at least much larger than il^, as we shall quantify in 
Sec. IV A) , in which case the dynamics of the cavity mode 
can be ignored or adiabatically eliminated. Note that the 
following analysis is also valid for the dark port fields en- 
tering and leaving an equal-arm Michelson interferometer 
with movable end mirrors and the differential motion be- 
tween these mirrors — but then the mirror mass m in the 
following discussion has to be substituted by the effective 



mass m/2. Then the quantum measurement process can 
be described by the following spectral densities [iO] 



2. 

ozz = tan C H n 

SpF = oi^ + 2mh , 



2h 
m fll 



SzF = fi tan C , 
where the coupling constant is defined by 



(48) 

(49) 
(50) 



(51) 



Here P is the circulating laser power and loq the laser 
angular frequency. 

In the free- mass limit, i.e. with ujm 0, the condi- 
tional variances simplify to 



\/2mn. 



1 + tan2 C + X 



(l + 2^2)(l + tan2C + 2e2)-tanC 



x/2 



Vxp 



1 + 2^2 )(1 + tan2 C -I- 2^1) - tanC 



(l + 2C2)(l+tan2C + 2e2)_tanC 



(52) 



(53) 



(54) 



while the purity of the state is given by 



det V = - ((1 + 2e|) (1 + 2^2) + tan2 C) 

>^ {l + 2^F^,f . (55) 

Here we have defined the two ratios, = ^p/^q and 
S,x = ^q/^x with riq = a/^/hm. Then we can recover 
Eq. (HI) and with the AppendixiBl wc obtain 



A/'cff 



2 5'sQL(f^c 



(56) 



In the quantum-noise limit with £,f = £,x = Q the con- 
ditional state is pure for any measurement frequency fJg 
and homodyne detection angle C. This defines the con- 
ditional ground state of a free mass 



^xx 



rig 



V2 



V 

Vxp - 2 



(57) 
(58) 
(59) 
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Here the Wiener filter functions for position and momen- 
tum become equal to simple decaying cosine functions at 
the measurement frequency and read 

K,{t) = V2nq e--^ cos (^^^ , (60) 

Kp{t) ^ Vlmnle-'-^ cos + (^1) 

Then it becomes apparent that the inverse of the mea- 
surement frequency 17, (at which the total noise ap- 
proaches SQL the most) is also the time scale at which 
information regarding test-mass position and momentum 
have to be extracted from the output data. 

It is clear from Eq. ([55)1 that one should measure the 
phase quadrature, i.e. C = 0, in order to minimize the 
uncertainty product. That means that in order to obtain 
a small uncertainty product, it is not required to remove 
the quantum back-action noise from the output, as one 
would benefit from when trying to detect gravitational 
waves — to the contrary, for that purpose it would be 
even destructive to do so. This is understandable, since 
here the aim is to learn as much as possible about the 
mirror motion, and the effect of quantum back action is 
an important content of the mirror motion. For example, 
suppose we chose to measure a back-action-evading ob- 
servable using an oscillator with negligible classical noise. 
In this case, we would have an output channel that has a 
minimal power spectrum, and is hence ideal for measur- 
ing any classical force that acts onto the mirror. Fluc- 
tuation in the mirror motion and the momentum driven 
by back action, on the other hand, would almost diverge 
around the resonant frequency. The output field, con- 
taining absolutely no information about the back action, 
would not be able to remove this fluctuation via condi- 
tioning. This would then result in a conditional state 
equal to the unconditional state with very large vari- 
ances. 

We learn from Eqs. ((5^ - ([5i)) that the effect of the clas- 
sical force noise in the conditional variances is suppressed 
with a higher measurement frequency (random force has 
less time to act and accumulate) while the classical sens- 
ing noise is suppressed with a lower measurement fre- 
quency (random sensing noise has longer time to average 
out). Moreover, in the absence of any classical sensing 
noise the test-mass state becomes pure with an infinitely 
strong measurement (in which Qq approaches -l-oo): all 
classical forces acting on the test mass can be neglected 
in presence of the strong back-action force and the test 
mass reaches the conditional ground state at the infinite 
measurement frequency. Vice versa, in the theoretical 
absence of classical force noise, the test-mass state be- 
comes pure in the limit of an infinitely weak measurement 
(in which fi^ approaches 0): if the test-mass motion is 
only driven by the measurement's back action but this 
motion is then unfortunately hidden in the measurement 
output because it is covered by the classical sensing noise, 
the best idea would be not to measure the test mass at 



all. If both classical noise sources are present, the un- 
certainty product is minimized further with an optimal 
power which accomplishes a balancing between classical 
force and sensing noise, i.e. with = S,x- This pro- 
duces an equal sign in the second line of Eq. (|55[) and 
is true for a measurement frequency of fig = ilci- This 
simply means that the quantum noise should touch the 
free-mass SQL at the frequency where the classical noise 
has the maximal separation to that limit (cf. Fig.[2|). 
The expression of the minimal uncertainty product and 
the one of the minimal effective occupation number are 
then functions of the classical-noise-SQL-beating factor. 

For a mechanical object there is not a fundamental 
definition of which state is vacuum, since it is possible to 
vary the potential well it lies in. Nevertheless, states with 
non-zero correlation among displacement and momentum 
can always be regarded as squeezed, and we can often 
discuss whether one state is more position squeezed (mo- 
mentum anti-squeezed) or more position anti-squeezed 
(momentum squeezed) compared to another. We have 
illustrated the squeezing situation, by plotting noise el- 
lipses obtained with various homodyne detections (see 
upper panel of Fig. [3]). As the homodyne detection angle 
C vary from optimal value to — 7r/2, the semimajor axis 
of the noise ellipse, i.e. the anti-squeezed quadrature, be- 
comes rotated into the direction of the position. At the 
same time the ellipse becomes more stretched, i.e. the 
squeezed quadrature is more squeezed, while the anti- 
squeezed quadrature is more anti-squeezed. As C vary 
from to 7r/2, the semimajor axis of the noise ellipse 
also rotates into the direction of the position but the el- 
lipse becomes rather bulged. Furthermore, we show how 
the noise ellipse changes, as the measurement frequency 
is not chosen to be equal to the geometrical mean of the 
two classical noise frequencies, i.e. fi^ 7^ f^ci- For phase 
quadrature detection as in the lower panel of Fig. [31 a 
slow measurement, i.e. with a low measurement fre- 
quency flq < ricii generates a position anti-squeezed (and 
momentum squeezed) conditional state; while a fast mea- 
surement with flq > fici generates a position squeezed 
(and momentum anti-squeezed) conditional state. In ad- 
dition, the deviation from the optimal measurement fre- 
quency, i.e. Clq 7^ Qci, always produces a less pure states. 



C. Very low finesse cavity and free mass scenario 
witii squeezed vacuum input 

Up to now we have only treated in-going coherent vac- 
uum states but one could also think about squeezed vac- 
uum states coupling to the mirror Q. This corresponds 
to inserting squeezed states into an interferometer's dark 
port. By doing so the quantum limited sensitivity of an 
interferometer can be enhanced 1491 . The free-mass con- 
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FIG. 3: Test-mass squeezing normalized with respect to the 
conditional ground state of a free mass for optimal measure- 
ment frequency Q,q — Slci = VTlpTl^ and different homodyne 
detection angles (upper panel) and for different measurement 
frequencies at phase quadrature detection (lower panel) in- 
cluding in both cases a certain classical noise budget: we have 
chosen fl^/ilF = 5. 



ditional state is given by 

h 



14: 



(A^ + 2^2 )(A?. + 2^2) _ sin2^op sinh2rop 



(62) 



(A^ + 2^j,){\i + 2^2) „ sin 2^op sinh 2ropy , 

(63) 



- I ( \/(^+ + ^^'f)(^- + 2Cx) - sin2^op sinh2rop 



cosh 2rop ±cos 2(pop sinh 2rop, 



Here we have defined A^ 
where (20/ In 10) Top > gives the optical squeezing 
strength in dB at a squeezing angle of iy9op- Then the 
purity of the conditional state can be inferred from 



det V = — ((A^ + 2e|)(A^ + 2^1) - sin^ 2^op sinh^ 2rop) 



>^[i+^eFe, + 2{eF+ cosh2rop 

- 2 - ej sinh2rop) > ^(1 + 2^^^.)' 



(65) 



The equality of the first inequality sign in Eq. (j65p is 
achieved at ipop = for > Q and at ipop = 7r/2 for 
< Cxj i-S- by squeezing either the phase or the ampli- 
tude quadrature, respectively. Note that in Eq. ((65| for 
any Jig, the same minimum as in Eq. (j55p is reached if 
Top = arctanh(|5|, — (^\/ {(^p + ^2))/2 — even when hav- 
ing S^p ^ £,x- Therefore, even with input squeezing, the 
conditional state cannot become more pure than with 
coherent input, but the demands on the required mea- 
surement frequency — and with this the constraints on 
the optical power, which is needed in order to obtain a 
certain uncertainty product, can be relaxed. In real ex- 
periments the optical power is of course always limited 
and squeezed input becomes an very important tool. 

It has turned out that the conditional variances are 
in principle even analytically equivalent in the follow- 
ing two cases: (i) input-squeezing at a flexible but 
frequency- independent angle or (ii) flexible amount of 
available optical power and a flexible but frequency- 
independent homodyne detection angle. This can eas- 
ily be seen by replacing the homodyne detection an- 
gle tan C — > sin 2(y9or> sinh 2ron and the measurement fre- 



(64) 



^op Olllll ^1 op 

quency Vlq — + A+fi^ in Eqs. 152]) - (|54)) . Then we simply 
end up with Eqs. Here we can directly see that 

using input squeezing allows to change the parameters 
within such as the optical power, the laser frequency 
and the mirror mass but by modifying the input squeez- 
ing parameter A+ we can at the same time maintain the 
measurement frequency. 

Even though a homodyne detection different from the 
phase quadrature and input squeezing do both not help 
with increasing the purity of the conditional state they 
increase the squeezing of the conditional test-mass state. 
Furthermore, with a certain homodyne detection angle 
or with a certain input squeezing it is possible to mini- 
mize the position and momentum correlation in the con- 
ditional state. 



IV. MACROSCOPIC ENTANGLEMENT 

The concept of entangled states is one of the most 
important phenomenons when entering the quantum 
regime. In Sec. lV Al we will see that the entanglement 
between the cavity mode and the mirror motion is in fact 
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responsible for a degeneratfon of the purity of the mirror's 
quantum state. In this section we will explain how we can 
prepare quantum entanglement in position and momen- 
tum between the centers of mass of the two end mirrors in 
the north i" and the east x° arm of a simple, but power- 
recycled Michelson interferometer using the conditional 
states as derived in the previous section. The end mirrors 
are suspended as pendulums but with a very low eigen- 
frequency. Therefore, such an experiment would be in 
direct analo gy t o the Einstein-Podolsky-Rosen gedanken 
experiment [50| . Here we will basically extend the dis- 
cussion carried out in Ref. 14011. 



with 



V V — 

* nil * 



(V + yd )/2 V +V'^ 

\ \ * xp ' ^xpJi PP pp j 

^y^XX ^Xx) I ^ ^K^Xp ^Xp) 1^ 
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Recall that the common (x" -I- a;") and the differential 
— x'^) mode of motion between the two end mirrors 
are independent and can each be sensed by a homodyne 
detection at the bright and the dark port, respectively, 
as suggested in Ref. [i^. As already mentioned before, 
if using the reduced mass of the mirrors, the conditional 
variances as derived in Sec. lIIII hold for the differential 
mode observed at the dark port of the interferometer. 
But in order to describe the common mode they have 
to be slightly modified since the power-recycling cavity 
— here with high bandwidth and therefore an adiabat- 
ically eliminated cavity mode — with transmissivity r 
enhances the measurement strength o!^ = 1It a> a and 
is therefore different to the one associated with the dif- 
ferential mode OL^ = a. Furthermore, the common mode 
will suffer additionally to the classical force noise and 
the classical sensing noise — note that we suppose that 
these two classical noise sources are equally distributed 
into common and differential mode — from laser noise, 
since the in-going modulation fields at the bright port are 
usually not in vacuum states. We have to make the fol- 
lowing additional replacements in Eqs. ([5^ - (|M|) in order 
to obtain the conditional variances of the common mode: 
(1 + 2^2) ^ ^Sa,a, + 2^2) and (1 + tan^C + 2^1) ^ 
(S'ajaa + 5*0101 tau^ ( + 2^^) as well as replacing the de- 
tached (— tanC) in each variance by {—Saiai tan^). Here 
S'aiai, Sa2a2 ^ 1 ^-rc thc (frcqucncy independent) spec- 
tra of the technical laser noise in amplitude and phase, 
respectively. 



Then we can assemble the conditional state of 
the entire system: the combined covariance among 
simply reads 



V. 



total 



V„e V„ii 



(66) 



This combined covariance matrix is very similar to the co- 
variance matrix for the amplitude and the phase quadra- 
ture of two output light beams which have been created 
by overlapping two continuous Gaussian light beams on 
a beam splitter. Note that overlapping two light beams 
which are differently squeezed in amplitude and phase 
quadrature on a beam splitter is a very common way of 
how continuous variable entanglement is created in op- 
tics [m, [13 • In the mirror case the common and the 
differential mode are mathematically overlapped to give 
the motion of each individual end mirror. 

In the following we will use the logarithmic negativity^ 
which was introduced for an arbitrary bipartite system in 
Ref. [s^, as a quantitative measure of the entanglement. 
For our state it reads 



£V = max[0, — logj 2<t /ft] 



(67) 



VS^ -4det V)/2 and S 



yc yd _|_ yc yd 
XX pp ' pp XX 



- 2^ yd _ 

xp xp 



where we have a = y (E 
det ynii det V^c - 2 dot Kc 

The higher the value of Ej\/ the stronger the entangle- 
ment. 

Recall that there exists a frequency band with sub- 
SQL classical noise if Vlx/VLp > 2. However, thc ex- 
istence of entanglement sets a slightly higher threshold 
value for this frequency ratio depending on the strength 
of laser noise as it is shown in Ref. [40| . We know from 
Sec. lIIII that the uncertainty product of each individual 
mode — common and differential — is minimal for a 
phase quadrature detection ('^•'^ = but it has turned 
out that this is not the optimal choice for the prepara- 
tion of entanglement. If the homodyne detection angle 
approaches — 7r/2, each mode can become more squeezed 
depending on the measurement frequency. 

The entanglement between the two mirrors — created 
by overlapping two modes — increases with the squeezing 
of the individual modes and with the angle separating the 
squeezed quadrature of the two modes. Then it is obvi- 
ous that one should not observe common and differential 
mode via phase quadrature detection, but that there is 
a certain value, for each the common and the differential 
mode, of fi^'d and —tt/2 < C,^''^ < which is optimal for 
the entanglement and maximizes the logarithmic negativ- 
ity (cf. the solid line in Fig. [3]). These optimal parameters 
depend of course on the classical noise, but are usually 
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FIG. 4: Logarithmic negativity versus Q.^IQ,f maximized 
with respect to f2q and using phase quadrature detection 
(dashed red line) as well as additionally maximized with re- 
spect to C,'^ and C,^ (solid purple line). In both cases no laser 
noise is assumed, i.e. with Sa^ax = Sa-^a-^ ~ 1- At some po- 
sitions optimal parameter values for Q'^, ^'^ and are 
given in the plot. 



characterized by a high measurement frequency and a 
detection close to ~tt/2 for both modes. That means 
that the states are totally driven by radiation pressure 
which is in turn monitored by reading out close to the 
amplitude quadrature. 

Furthermore, we have found that the laser noise — 
entering at the bright port and in a totally balanced in- 
terferometer only affecting the common mode — can the- 
oretically be almost suppressed with the optimally high 
measurement frequency and the optimal homodyne de- 
tection angle. Therefore, the resulting maximal entangle- 
ment — represented by the solid line in Fig. 5] — is inde- 
pendent of the laser noise. But the parameters — such as 
optical power and the fine-tuning of the homodyne detec- 
tion angle, which are required in order to reach the max- 
imal entanglement, are different for different strength of 
laser noise and are far away from any realistic situation 
(cf. the dots in Fig.|l]). 



V. INTERFEROMETERS WITH 
NON-MARKOVIAN NOISE 

A. Cavity with finite bandwidth 

In realistic experimental situations, the noise sources 
are usually not Markovian. In this section we want to 
generally study quantum state preparation in the back- 
ground of non-Markovian noise sources. We start with 
allowing the quantum noise spectral density to become 
frequency dependent. 

If we consider a cavity of length L with a finite half 
cavity bandwidth 7 and a movable end mirror, the quan- 
tum noise is indeed non-Markovian and the Hcisenbcrg 



equations of motion in frequency domain modify to 

in 



y 



sinC 



7- 



7-if7 

J -in 



ai + cos C X 



02 



in h 



m{n'^ + Z7„iO - uj^) 



"f-in 



aai + 



(68) 

m 



where we have approximated exponential functions by 
rational functions of n. Now the measurement frequency 
becomes 



'^1 mhL"f 



(70) 



which corresponds to the frequency where the quan- 
tum noise spectral density of the associated adiabatically 
eliminated system touches the SQL. If the phase quadra- 
ture of the outgoing light {( = 0) is detected, the eight 
zeros of its spectral density Syy are given by ±ai ± i6i 
and ±02 ± 162 with (for simplicity for a free mass, i.e. 

7m = UJ.,n = 0): 



(71) 
(72) 




where 



1 + 1. The zeros are 



required for the spectral decomposition introduced in 
Secini After straightforward algebraic manipulations, 
one arrives at the following conditional second-order mo- 
ments: 



h'-f 



PP 



6m(f]^^'^)2 

"120(17^^ 
fi7^ 



{cl + 3cl + 3ci + C3) , 



(73) 



(3c5-|-15c;[ + 20cJ + 60c3 + 60c5), (74) 



c?(ci + 2)2, 



where the coefficients are given by 
2 

c„ = -5 [(ai + i5i)" + (02 + 162)" 



(75) 



(76) 



The conditional variances given in Egs. ([75)1 - ([75]) . al- 
though slightly complicated, are still analytic; we can 
still draw some important conclusions from these second- 
order moments. By expanding the quantity U from 
Eq. ([31) in terms of n"^^^ /j, we obtain 



U ^l- 



1 n 



2V2 1 



^+0((f]-77)' 



(77) 
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FIG. 5: Test-mass uncertainty product U (solid line) com- 
pared to the entanglement between test mass and cavity 
mode (dashed line) both versus the dimensionless ratio £7^^^ / 7 
through which both quantities are totally described. No clas- 
sical noise is present. Free mass limit is used, i.e. 7™ = uom = 
0, and the phase quadrature <^ = is detected. Entanglement 
is quantified by the logarithmic negativity. 




FIG. 6: Contour plot of the normalized test-mass uncertainty 
product given by f/ = 2/hyJVxxVpp — V^p versus Q.q'^^ /j and 
•jJmh for quantum noise only. Again phase quadrature (^ = 
0) is detected. 



which reveals that, even in the quantum noise limited 
case, the conditional state of the test mass cannot be 
pure as long as > 0. This is in contrast to the 

Markovian Umit (7 — * 00) discussed in Scc. lIIIl where 
the conditional state is always pure in the absence of 
classical noise [cf. Eg. ([551) ]. Fig. O shows the purity of 
the test mass versus the dimensionless ratio 57^^^ /7. 

In the case of a finite cavity bandwidth, the light is 
stored inside the cavity for some time. The information 
carried by the light concerning the test masss state can- 
not leave the cavity instantaneously and hence is not ac- 
cessible for the conditioning process. Consequently, the 
intra-cavity mode needs to be taken into account for a 
complete characterization of the system. The residual 
second-order moments required for completing the corre- 
sponding (4 X 4) conditional covariancc matrix can be 
obtained in the same way as Egs. ([75)1 - fTS]) . It turns 
out that the composite system is indeed a pure one even 
though each individual system resides in a mixed state. 
This is a clear evidence of entanglement between the con- 
ditional states of the test mass and the cavity mode. 
We have also plotted the logarithmic negativity in Fig. (5] 
The test-mass state's purity decreases and test-mass- light 
entanglement increases with smaller bandwidth and with 
higher measurement frequency fi^^'^. Note that the un- 
certainty product as well as the logarithmic negativity 
do not diverge. But Fig. [5] indicates that, as long as 
riq*^^ <SC 7, we can neglect this effect and adiabatically 
eliminate the cavity mode as performed for the power- 
recycling cavity in Scc. lIVI 

Fig. [5] further shows, that the purity increases with 
higher mechanical eigenfrequency tj„i, depending on the 
measurement frequency. Let us consult the following 
hand- waving argument: with increasing uj^a the mechan- 
ical oscillator and the optical oscillator, which would res- 
onate at modulation-frequency zero, become more sepa- 
rated in the frequency space and therefore their entangle- 



ment decreases. And with decreasing entanglement the 
test-mass state becomes more pure. But the regime with 
such high mechanical resonance frequencies is usually not 
available in actual GW detectors. 

When including the two classical noise sources from 
our simple model the conditional state becomes more 
and more mixed with an increasing classical noise level 
as shown in Fig.[7l Note that this is of course not an 
effect of entanglement, the classical noise rather destroys 
the quantum entanglement between the test mass and 
the cavity mode. The purity also depends on the ra- 
tio between the measurement frequency and the optical 
bandwidth. Recall that the classical force noise increases 
with higher p while the classical sensing noise increases 
with lower O^,. If we only take sensing noise into ac- 
count we know for sure that the motion of the test mass 
is solely driven by quantum back-action noise. A high 
sensing noise level randomizes the measurement record 
and hence the delimitated accessibility of the intra-cavity 
field due to the finite bandwidth is insignificant. Conse- 
quently all curves roughly coincide for fi^^^ > fij, in the 
upper panel of Fig. [71 

Note that in the planned Advanced LIGO detector 0] , 
the measurement frequency $7^'^'^ is planned to roughly 
coincide with the half cavity bandwidth at 7/(27r) ~ 
100 Hz — corresponding to the solid lines in Fig. [71 Fur- 
thermore, we expect the suspension thermal noise to have 
a 51i?/(27r) ~ 30 — 40 Hz — that would be less than 
j-jcav^2 — but the coating thermal noise may provide 
a ria: that only coincide with or is just marginally 
higher. We can infer form Fig. [71 that the quantum state 
of the interferometer's differential mode is mainly con- 
strained by classical sensing noise which entails a lower 
bound of f/ > 5 and gives an effective occupation number 
of A/'cff ^ 2. For a more detailed discussion see Sec. lVIl 



13 



: — n-Vr = 


0.1 // /' 




1 // / 




10 // / 










'''■^'^with classical (white) sensing noise 



0.1 0.2 



0.5 



1 



10 



'^^Vy = 0.1 


/ / / 


— n-Vy = 1 


* / // 


... n^^Vy = iO 


/ V 


— n^^Vy-0 / 






/ .// 












with classical (white) force noise 



0.1 0.2 



0.5 



1 



10 



FIG. 7: Test-mass uncertainty product I] versus classical 
sensing noise and without force noise (upper panel) as well 
as versus classical force noise without classical sensing noise 
(lower panel), both for different examples of ratios between 
measurement frequency and bandwidth. Free mass limit is 
used, i.e. a;™ = 0, as well as phase quadrature detection 

C = o. 



B. Detuned cavity 

A cavity which is detuned by A from the carrier's fre- 
quency makes the power inside the cavity also dependent 
on the motion of the test-mass mirrors. This creates an 
optical spring [s^ or an optical anti-spring both shift- 
ing the (free) mechanical and the (free) optical resonance 
frequency in the complex plane. In future GW detectors 
such as Advanced LIGO Q , the optical spring effect will 
be used to up-shift the real part of the mechanical res- 
onance frequency into the detection band. Recall that 
the optical spring as well as the optical anti-spring usu- 
ally introduces instability to the system which has to be 
cured with an appropriate linear feedback control [i^ . 
But it is straightforward to show that the conditional co- 
variance matrix does not change under any ideal, linear 
feedback control. 

The Heisenberg equations of motion for such a system 
can be found in many previous works, see e.g. Eqs. (39)- 
(41) of Ref. [5^. Unfortunately, analytic expressions 
for the conditional covariance matrix are cumbersome, 
and we only report numerical results. We also restrict 
ourselves to quantum noise, and reading out the phase 
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FIG. 8: Contour plot of the normalized test-mass uncertainty 
product given \sy U = 2/h^VxxVpj, — Vip versus the ratio 
between detuning A and bandwidth and Qq'^^ for quantum 
noise only. Free mass limit is used, i.e. uJm = 0, as well as 
phase quadrature detection ( = 0. 



quadrature, i.e. C 



that the measurement frequency fi^^^ 



(cf. Eq. (3.2) of Ref. \4$). Note 
is defined in the 

same way as for the tuned finite bandwidth configura- 
tion (cf. Eq. dZni)). 

Fig.[S] shows that detuning a cavity from the carrier 
frequency properly, can increase the purity which comes 
from the fact that the quantum entanglement between 
test mass and cavity mode is decreased. In the regime 
of a blue detuned cavity (A > 0) — producing an opti- 
cal spring — and for fi^^'^ < 7, Fig. [8] simply agrees with 
Fig. [HI Here at fixed measurement frequency fi^^^ < 7 
a higher detuning A gives a less shifted mechanical res- 
onance, so-called optomechanical resonance, and at the 
same time it corresponds to a higher optical resonance. 
Therefore, again the two oscillators are more separated 
in the frequency space and their entanglement decreases. 
Interestingly, for higher fl"^^"^ the test-mass state could lo- 
cally appear more pure in the red detuned cavity regime, 
i.e. at a certain A < 0, which produces an optical anti- 
spring. Note that the uncertainty product diverges for 
an infinitely red detuned cavity (A — > —00). For these 
facts we unfortunately have not found any intuitive ex- 
planation. This will be a subject of further investigation. 



C. Non-Markovian classical noise 

In the following we consider a more realistic example 
configuration involving multiple colored classical noise 
sources. The classical noise contributions are highly non- 
Markovian and they tend to rise fast in the low frequency 
regime, which is ignored by a simple Markovian noise 
model. First we restrict ourselves to an idealized noise 
budget of an advanced interferometric gravitational wave 



14 




f [Hz] 

FIG. 9: Upper panel: seismic noise pre-estimated by the 
simulation tool Bench [s^ (red) and a fit by rational func- 
tion (gray). Lower panel: straw man (colored) classical noise 
budget of an advanced interferometric gravitational wave de- 
tector: Seismic (violet), suspension thermal (blue, follows 
~ 1//^/^) and internal thermal (green, follows ~ 1//^''^) 
noise spectra are shown as well as the total noise (black) — 
including Markovian quantum noise with phase quadrature 
readout. For the thermal noise sources we employed a Fade 
approximation that generated the 1// spectra between 0.1 Hz 
and 2 kHz. 



detector, shown in the lower panel of Fig. [^1 Here only the 
dominating force and sensing noise sources are consid- 
ered. Additionally we assume that the gravity gradient 
noise can be suppressed completely through monitoring 
the ground's motion. Especially seismic noise dominates 
the entire spectrum below 10 Hz. In order to apply the 
numerical Wiener filter procedure (cf. Sec.|TT|, all clas- 
sical noise spectra need to be approximated by rational 
functions of O^. This is illustrated by Fig. [^K upper panel) 
where the seismic noise spectral density, pre-estimated 
by the simulation tool Bench [s^ . is approximated ac- 
cordingly. The seismic noise spectrum is constant below 
0.25 Hz, then it drops as ^ 1//^ between 0.25 Hz and 
2 Hz and finally it drops as ^ 1//^° above 2 Hz. The 
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FIG. 10: Uncertainty product of the conditional state versus 
sensing noise cut-off frequency for three different sensing noise 
levels. Force noise contributions are the same as in Fig.|9] 
The second-order moments formally diverge when downshift- 
ing the cut-off frequency. 



suspension thermal noise constitutes a second force noise 
contribution which drops as ~ 1/ /^^^ above the pendu- 
lum eigenfrequency (at 1 Hz) and it intersects the SQL at 
20 Hz. Such a frequency dependance presumes structural 
damping. Above 3 Hz the internal thermal noise follows 
^ 1//^/^ and it intersects the SQL at 500 Hz. We have 
employed the Fade expansion in order to simulate the 
behavior of the spectral densities. 

It should be emphasized that the conditional second- 
order moments can diverge, if the sensing noise rises to- 
wards low frequencies, and therefore a cut-off frequency 
must be chosen carefully. This issue is illustrated by 
Figs.fTOl where the cut-off frequency of the sensing noise 
is varied, while the classical force noise contributions are 
held fixed. This divergence can be explained as follows: 
for a free mass, the effect of radiation pressure noise di- 
verges towards low frequencies. Hence the boundedness 
of the conditional variances depends crucially on the the 
motion of the test mass at low frequencies. Furthermore 
the mirror thermal noise with S{f) ~ 1// formally leads 
to a logarithmic divergence. In the case of a real exper- 
imental set-up the mirrors are suspended as pendulums 
and the mirror thermal noise should also exhibit a low- 
frequency cut-off — but more importantly the low fre- 
quency noise will be canceled out in a subsequent verifica- 
tion stage, as it will be shown in a forthcoming paper [43 | . 
Such a cancelation arises from the fact that for the mirror 
thermal noise with frequencies lower than the inverse of 
the sum of the preparation-stage and verification-stage 
measurement time scales, their contribution to errors in 
the preparation and verification measurements are the 
same, and therefore cancels out, when the two sets of 
data are compared with each other. This argument jus- 
tifies an increase of the cut-off frequency to a level of 
around 3 Hz. 
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VI. ADVANCED LIGO CONFIGURATIONS 



In this section wc will investigate the performance of 
the planned Advanced LIGO detector Q , a second gener- 
ation gravitational-wave observatory, towards the prepa- 
ration of test-mass quantum states. It is planned that 
this large-scale laser interferometer (cf. Fig.[T]) — 4 km 
long arm cavities consisting of 40 kg mirrors — starts its 
operation in 2014. It will be nearly quantum noise lim- 
ited in most of its frequency band (10 Hz to 10 kHz), and 
will operate near or at its SQL. Our previous investiga- 
tions in Sec. mil and IVl have suggested that such a SQL- 
sensitivity allows to prepare nearly Heisenberg-limited 
quantum states of macroscopic test masses. Note that 
we will consider the differential mode of the interferome- 
ter's four movable arm-cavity mirrors which is equivalent 
to a single movable mirror in a single detuned cavity (SGj 

— with one quarter the mass of each individual mirror, 
i.e. 10 kg. 

The classical noise budget of the Advanced LIGO 
detector has been estimated by the simulation tool 
Bench [11]. We choose the same type of spectra as for 
the example configuration in Sec. IV C[ i.e. with iden- 
tical power laws and cut-offs in frequency, and adjust 
the parameters such that the predicted Advanced LIGO 
classical noise budget is well approximated. In contrast 
to Sec. IV CI Advanced LIGO comprises finite bandwidth 
cavities, which gives rise to a non-Markovian quantum 
noise. Moreover, the detuned signal-recycling technique 
introduces even the optical spring into the dynamics of 
the mirrors just as in the case of a detuned cavity (cf. 
Sec. IV Ap . As an example, the quantum noise of the Ad- 
vanced LIGO broadband configuration, which is optimal 
for the detection of neuron star binary inspirals, is plot- 
ted in the upper panel of Fig.fTTl 

We have carried out a full parameter search over the 
space of signal-recycling parameters in order to optimize 
the configuration with respect to the uncertainty product 
of the conditional state. We have fixed the characteris- 
tic frequency tc (as defined in Eq. (20) of Ref. [5^) of 
the system, which is basically determined by the fixed 
circulating optical power of 800 kW which in turn de- 
termines the measurement frequency. Furthermore, we 
have only considered a homodyne detection at the Ad- 
vanced LIGO broadband configuration quadrature, i.e. 
C = 0.7 TT. The result of this optimization is shown in the 
lower panel of Fig.[TT] which depicts the effective occu- 
pation number ACff — as introduced in the AppendixiBl 

— versus the effective detuning A and the effective band- 
width e — these two quantities arc defined in Eq. (18) of 
Ref. [5^. It clearly shows that the purity of the condi- 
tional state benefits from a restoring optical spring, i.e. 
a positive detuning facilitates the preparation of macro- 
scopic quantum states as it was shown before. Note that 
increasing the effective bandwidth e always gives rise to 
an additional improvement, which has also been clarified 
before. The Advanced LIGO broadband configuration 
with A = 27r290Hz and e = 27rl20Hz is marked with a 
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FIG. 11: Upper panel: spectral densities of main noise sources 
present in the Advanced LIGO detector: seismic, suspension 
thermal and internal thermal noise (gray lines) as well as two 
examples of the quantum noise (magenta and purple lines). 
Pre-estimated classical non-Markovian noise budget is fitted 
by rational functions with characteristic spectra as in Fig.|9] 
Lower panel: effective occupation number A/'efi of conditional 
state (of the differential mirror mode) versus effective detun- 
ing A and effective bandwidth e. The dots mark the Advanced 
LIGO broadband configuration state (magenta) and lowest 
occupation number state (purple) when ^ — 0.7 tt. 



dot in the lower panel of Fig.fTTland gives AfeS « 2.2 while 
the other dot marks the purest state at A = 27r 500 Hz 
and e = 27r400Hz which gives JVcS « 1.9. An additional 
optimization of the homodyne detection angle decreases 
this number only marginal. 

Aside from the currently estimated classical noise bud- 
get, a more optimistic scenario [57| has been investigated, 
in which the seismic and suspension thermal noise are re- 
duced by a factor of ten, while the coating thermal noise 
is lowered by a factor of three (in amplitude). Here the 
cut-off frequencies remain the same. A rough optimiza- 
tion has revealed that the minimal achievable effective oc- 
cupation number drops down to Afcs ~ 0.38 (cf. Fig. fT^ 
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A/(2 n) [Hz] 

FIG. 12: Upper panel: spectral densities of main noise sources 
present in an improved Advanced LIGO detector. Lower 
panel: effective occupation number A/'eff of the differential 
mirror mode's conditional state versus effective detuning A 
and effective bandwidth e for phase quadrature detection. 
The dot marks the lowest occupation number state. 

for phase quadrature detection. Further major improve- 
ments regarding the classical noise level in the Advanced 
LIGO detector — especially the laser noise — will even 
allow to entangle the cavity mirrors in the north with 
those in the east arm. 



VII. CONCLUSION 

This paper has been devoted in great detail to a survey 
of the first principles in the preparation of macroscopic 
Gaussian quantum state of non-massless objects. We 
have motivated and introduced the Wiener filter method 
in this context — as an advantage over the stochastic 
master equations — and have given a simple analytical 
expression for the covariance matrix of a system under 
any continuous linear Markovian measurement process. 



We have shown that in absence of any additional noise, 
the conditional state is totally determined by the mea- 
surement noise. Moreover, the purity of the conditional 
state is even equal to the purity of the underlying mea- 
surement process (cf. Eq. (|4ip ). This provides an impor- 
tant insight into the understanding of conditional states 
which was probably not communicated before. 

In Markovian measurements with non-correlated shot 
and radiation-pressure noise, we have shown that the ef- 
fective occupation number of the conditional state is con- 
nected to the factor the device beats the SQL, i.e. to the 
bandwidth within which the classical noise is below the 
SQL (cf. Eq. gT]) and Eq. ([51])). For Sd < S'sql around 
the frequency Sid, at which the two classical noise spec- 
tra intersect, we find A/'eff < 1/2. We have shown that 
neither a balanced homodyne detection of a non-phase 
quadrature nor input-squeezing would help to get a more 
pure state — but they can significantly steer the shape 
of the conditional state, e.g. the test-mass squeezing. 

Furthermore, we have motivated that a simple power- 
recycled Michelson interferometer is the ideal device to 
prepare macroscopic entanglement [40j . We have shown 
that the existence of entanglement in position and mo- 
mentum between the two end mirrors is closely related 
to the factor at which the classical noise beats the SQL: 
a quantum measurement with a fiexible but frequency- 
independent homodyne detection angle and no restriction 
to the optical power as an example, theoretically requires 
the classical noise to be at least a factor of 1.5 below the 
free-mass SQL at a certain sideband frequency. 

Moreover, we studied mirror quantum-state prepara- 
tion in non-Mar kovian quantum- measurement systems. 
In the first instance we have considered the conditional 
quantum state of a test mass inside a finite-bandwidth 
system. It has been demonstrated that even a quantum 
noise limited configuration does not allow the prepara- 
tion of a minimum Heisenberg uncertainty state — due 
to quantum entanglement between the test mass and the 
cavity mode, which has a non-zero lifetime. 

It has been pointed out that the purity of a conditional 
quantum state of macroscopic test masses can benefit 
from introducing an optical spring. This has been verified 
numerically for the quantum noise limited regime. 

Furthermore we have optimized the effective occupa- 
tion number of the differential mode of the planned Ad- 
vanced LIGO gravitational- wave detector in the presence 
of prc-estimatcd realistic decoherence processes. It has 
been confirmed that already a moderately reduced clas- 
sical noise budget, such as for an improved Advanced 
LIGO detector, allows us to prepare a nearly pure quan- 
tum state of the mechanical mode under consideration. 

Finally, we have explored the effective occupation num- 
ber achievable by Advanced LIGO, on the differential 
mode of its four macroscopic test masses. It has been 
shown that an occupation number of « 2.2 is readily 
reachable by the baseline design, a moderate shift in op- 
tical parameters can achieve w 1.9, while a moderate en- 
hancement in classical noise budget could achieve w 0.38. 
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Third-generation gravitational-wave detectors, or pro- 
totype interferometers specifically designed for testing 
macroscopic quantum mechanics would be able to sur- 
pass this moderate enhancement of Advanced LIGO, and 
reach deep into the quantum regime. 
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APPENDIX A: QUANTUM WIENER FILTER 



with 



-i^U,=o J = .xr'W (A2) 



If we then write the projection operator as a path inte- 
gral. 



[fi(t')=y(t)|t'<t] « 



i?[fc]e'^-=^''*''^(*')[*(*')-^(*''],(A4) 



we have 



J cx / D[k]tr 



p giE, ai«i(t)+i/!^ dt'k{t')[y(t')-y{t')\ 

„ (A5) 

For a Gaussian state p and any linear observable x = x\ 
if tr[/5 x] =0, then we always have 



Here we directly evaluate the conditional generating 
functional involving the observables xi (t) , 



J = tr 

= tr 



-cond iJ2i aixi{t) 



P 'Pmt')=yit)\t'<t] e' 



5Di aixi(t) 



J 



tr [p e"] = e 



(Al) Using this property, we obtain 



(A6) 



J cx y" D[fc]e"([^' ^iMt)+H^ rft'fc(t')S(t')]')/2g-i/I^ dt'k{t')v{t') 
Suppose Eqs. (dH) and ^ hold, i.e.. 



(A7) 



xi{t)=Ri{t)+ dt'Ki{t^t')y{t'), {Riit)yit')) ^0,W <t, 



then 



"^aixiit)+ I dt'k{t')y{t') 

1 J —CO 



L I 



where we have defined 



'^aiam{Ri {t)Ryn{t)) 



Ira 



dt'dt"k{t')k{t"){y{t')v{t")). 



OQ J ~OQ 



kit')=kit')+J2c^iKiit-t') 



(A8) 



(A9) 
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Using k(t) as the new integration variable, J can be re-written as 



J oc exp 



^^aiaMWRmit)) exp if dt' Y,aiKi{t - t')y{t') 

Ira J L I 



(AlO) 



which justifies Eqs. (fTi]) and ((15] 



APPENDIX B: EFFECTIVE OCCUPATION 
NUMBER 

The uncertainty product of Gaussian states is a true 
measure of the purity and therefore a reasonable measure 
of the quantum-ness of the state. Trying to reconstruct, 
as commonly done, the number of quanta, the so-called 
occupation number, may not always be the most funda- 
mental figure of merit: squeezed states, for example, can 
have high occupation numbers, yet they should be consid- 
ered probably more quantum than vacuum states. More- 
over, the definition of an occupation number requires a 
well-defined, real-valued eigenfrequency, which does not 
always naturally exist. 

The uncertainty product can be converted back into 
an effective occupation number by using the relation 

2 AAeff = I y^l4, V„ - 1/4, - 1 . (Bl) 

If a state has no correlation in position and momentum 
this effective occupation number should be interpreted 
as follows: suppose that the variances in position and 
momentum are given and produced by a perfect harmonic 
oscillator in a quadratic potential having an arbitrary 
but real- valued eigenfrequency locS- Then the effective 
occupation number is obtained by minimizing the total 



energy divided by the energy of each quanta with respect 
to that eigenfrequency Woff ■ This strategy reads 

AAcff ^ mm <^ + ^1 -of' (^2) 

where the minimum is achieved at 

Thus, the effective occupation number is the minimal 
occupation number one could obtain when assuming to 
have a harmonic oscillator with no correlation in position 
and momentum and an effective eigenfrequency as given 
in Eq. 

For a state with correlation in position and momen- 
tum, the effective occupation number still gives the min- 
imal occupation but with respect to two other orthogo- 
nal quadratures which are not position and momentum. 
Moreover, the effective occupation number is an interest- 
ing quantity because it in fact determines the von Neu- 
mann entropy of a state [sst as given by 

S = (AAeff + 1) log(AAeff + 1) - AAeft logA^eff ■ (B4) 
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